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Field of the Invention 

The present invention relate to a process for. increasing the efficiency ofl a 
computer in finite alenxonrt simulations by efficient automatic construction ofl 
suitable basis functions for computation of approximate solutions, and to a 
5 computer for pcrformtfnig tWat process. 

Background ofl the Invention 

A plurality of technical and physicall phenomena can be described by partial 
differentia! equations. They include: among others problems from fluid mechanics 
3G (for example, flow around an ainfoil]', electromagnetic field theory (for example, 



f 

- 2 - 

electrical field behavior in a transistor} or elasticity tli'cory (for example, 
deformation of a car- body). Accurate Knowledge and description of such' processes 
are a central element in the construction and optimization of technical objects. Tb 
save time-conisumiing and cost-intensive experiments., theno is great interest in 
5 computer-aided simulations. Finite element processes (FE processes) haw baoome 
established and have been tihle topic of intense research for. a long time. This also 
applies to automata mesh 1 generation processes as a foundation' for. construction of 
suitable basis functions. 

10 Figure 1 illustrates in the left half the prior art! in the piraoess of FE simulation for 
a linear boundary, value problem as a typical model example. Proceeding from the 
data describing the gjeometry of the teohnical object to be simulated, first a system 
of basis functions is constructed which! on the one hand enables fldfilmjent of 
boundary conditions and on the other is suited for approximation of thJe unknown 

15 solution. Tlien, using these basis functions a system of) linear equations is set up 
by numerical integration methods. Finally, the coefficients of the unknown 
approximation ape determined as the solution of this system ofl linear equations. 

FIE methods on tihoir use axe the subject matter of a series of patents. For example, 
20 U.S. Patent No. 4, 81)9^ 161 discloses a system where FE approximations of a large 
class of different iail equations are automated. U.S. Patient No. 5,731,817 discloses 
a process for generation of ltoxahedral mashes forming the foundation for a FE 
simulation process. 

25 In most EE processes af practical relevance basis functions are used which are 
defined on a decomposition produced by) generating a mesh of the simulation 
region. Fligurc 2a* shows a selection of conventional elements) their dimension, 
degree, smoothness and panarmetors aire listed in Figure 2b. A survey, of meshing 
mctilio.ds of planar regions ran be found fbn example in K. Ho-Le, Finite Element 

30 Mesh Generation Mcdhods: A. naviicwi and classification. Com. Aiided Design 20 
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£1988), 27 - 38. Generating ai mesh for complicated three-dimensional regions is 
extremely) difficult using the current stia.ta oil knowledge, as shown by S. Owen, A 
survey of unstructured Mesh Generation Technology, Proceedings, 7th 
International Meshing Round Table, Sand'ia National Lab' (1998], 239 - 2571. The 
5 processes require extensive amounts of computer time and are to some extent not 
yet completely; automated. But recently, there has been a series of very innovative 
new approaches. For example, AI. Fuchs, Optimierte Delaunayj-TriaTTgufierungen 
zur Vcrnetzung getrimmter NIURJBS-Korpar, University of Stiuiitgart, 1999, 
simulates a force distribution in ordeii to achieve an optimum distribution of 
1Q tiiiangulatian points. In U.S. Patent No. 5, 729, 670 two- and three-dimensional 
meshes are produced by solving flow problems; this is an interesting reversal 1, of the 
conventional FE mechanism. In addition, many algorithms have been developed to 
improve individual aspects of: mesh generation prooesses. For example DE 
19621 434 Al and U.S. Patent No. 5, 774, 696 describe a process for elimination of 
15 intersections with: prescribed edges or boundary surfaces in Delaunay 
fcriangulations. 

Mcshlcss FE-mctl>ods to daite have not acquired any t importance for applications. 
RotW in tbfc Lagrange multrplier method (see for example J.H. Bramble, The 
20 Lagrange Multiplier Method fbr Dimchlet's Problem, Matfh. Comp. 

37 £19ftl), 1-11), and also in the penalty method (see for example P. Boohev. and 
M. Gunaburgcr, Finite Elemont Methods ofi Least! Squares Type, SIAM review 
40 (11998), 7i89 - 837), the treaitmiant of boundary conditions represents a major 
problem in the use erf simple, stable basis functions. 

25 

Tn many technical simulations l automatic mcshl generation is very complex and 
requires by far the laTgcst part of the computer time. Furthermore, thte 
approximation power of the conventionally used linear and multilinear basis 
functions is low. To achieve accurate results, a large number of basis functions 
30 must he used, and thus, a correspondingly large system of equations must be 
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solved. Higher order trial functions on triangulatro.ns generally likewise Wave an 
unfavorable ratio between trio attainable accuracy and the number of basis 
functions used. Finally, smooth basis functions cannot be easily defined on 
unstructured mesh'es. Very special constructions are necessary] already for 
5 continuously difTerentiab'le elements (see Figure 2a). 

Summary of the Invention 

Th!e abject of the present invention is to increase the efficiency of known FE 
methods and computers which carry out FE methods by efficient construction of 
basis functions with favorable properties. In particular the meshing of the 

Id simulation regiion will be completely, eliminated, optional boundary conditions are 
fulfilled, accurate solutions are obtained with relatively few coefficients, and the 
resulting systerro of equations will he solvable efficiently. In this way the 
disadvantages of the prior, art wall' be overcome, and thus, the accuracy and speed 
of the simulation of physical properties in the engineering and optimization of 

15 technical objects will! be improved. 

Some central! terms and notations which are used in the following description of 
the process as alaimed in the invention will -be explained first. 

2Q The simulation region Q is a bounded set of dimension d => 2 ox d = 3 on which the 
physical quantities to be studied will be approximated by means of FE methods. 
The boundary of the simulation region is denoted by F. A grid with grid width h 
is defined as a decomposition of a subset of: the plane or the space in grid cells Z k . 
Depending on the dimension d each grfd cell is a square or a cube with edge length 

25 h. More precisely, 7,% = kh + % k] d , where h belongs to a set of integer d-vectors. 
Tlic uniform tensor product B-spBflnes in d variables of degree n with grid width k 
arc denoted by i^'scc for example 0. de Door, A Practical Ouide to Splines, 
Springer, 1978. T.hcy arc functions which can be continuously differentiated! (n— Uj 
times and which on the grid cells agree with polynomials of degree n\ as shown in 
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Figures 4a and 4b. Figure 4!ai shows flic support of flic B-splmc o.f degree 
7i = 2, dimension d^2 and smoothness nv = 1. In Figure 4b, flic resulting tensor 
product] B-spline b% is shown. The support Q,k> i.e. the union of all grid cells, on 
which 1 Che B-spline b k ils not identical to zero, consists of. (n + Tj d grid cells; more 
5 precisely, Q k = kH^ [0, (n + l)/i] rf . In all figures, the B-spline bn is marled at the 
point kh t i.e. for example in the case d = 2 at the lower left corner of the support. 
Eon EE simulations only thlose B-splines are important which have support 
intersecting the simulation region H; th'ey are called relevant B-splines. The 
relevant B-splines are again divided into two groupsj those B-splines for which the 
10 part of the support inside the simulation region is larger than a prescribed bound $ 
are called inner B-splioes. All other relevant B-splinies are called outer B-splines. 

Tllic abject of the present invention is achieved by tb!e process defined in claim L 
Special? embodiments of the invention are defined in the dependent claims. Claim 
15 11 defines a computer system as claimed in the invention. 

The right half of Figure 1 shows the incorporation of tine process of thle present 
invention into tiho course of a FE simulation i)n the prior art and. thle substitution 
of certain process steps of a EE simulation in the prior art by the process of thle 
20 present invention. 

Input 1 of the simulation region Q can be done via input devices, in particular also 
by storage of data derived from comipiiter-aidod engineering (OAD/|CAM)j. For 
example, the data used' in the engineering of a motor vehicle cart be incorporated 
25 directly into the EE simuilati'on 2n the present invention. 

In input 2 send storage nil the tiype of boundary, conditions natural and 1 essential 
boundary conditions arc distinguished. THe basis in the present invention 3s 
constructed for homogeneous boundary conditions of the same type. Tn particular., 
3d for essential boundary conditions, tifoo basis functions vanish on the boundary R 
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Infoomogeneous boiindlanyi conditions can he treated in the assembly of the FIE 
systems using methods which correspond! to the prior art. 

Finally, the control parameters are read ?n 3. They; rel'aco to t-We degree n andl title 
5 grid width 1 h of the B-splines to be used and the bound s for classifioationi of the 
inner- and outer B-splines. If specifications are omitted, all these input parameters 
oan be automatically determined by evaluation of merit functions which arc 
constructed empirically or analytically. 

10 Tha following construction of the basis functions of the invention is divided into 
the steps shown schematically in Figure 3, which wifl now be described. 

After reading in the simulation region fif, in the first process step a grid covering 
tike simnilation region Q is generated. Then it is checked which of the grid oells lie 

15 entirely iroside, partJallyi inside or noil inside the simulation region Q. The cell 
types 4 are determined, and this information about the cell types is stored. Thi9 
essentially requitoas inside/outside tests and determinations of intersections 
between the boundary; II of the simulation region Q and the segments or squares 
wJiicli hound the grfld cells. Figure 7 shows the input and output data for this 

20 process step. 

Tn the second process stop, using the information about the cell types, the relevant 
B-Bplincs aTc first determined Then the classification 5 into outer B-splines is 
performed; the corresponding lists of indices are denoted by / and ,K To this end, 
25 the size of tl>osc parts of the supports dfl the B-splines which lie wMiin the 

simulation region is determined' using the data obtained in the first prooess step, 
and compared with/ the prescribed Hound a. Fligure 9 shows the input and output 
data for this process step. 

30 In the third process step, coupling coefficients are computed 6; they join th'e 
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inner and outan B-spllines according to the rule 

Hence, an extended R-spiine is assigned to ea*;h inimer R-spline Bj. The 
construction and the properties of the index sets J{i) and the coupling ooefficienits 
5 e{j are given as follows. Tibe index sets ,7l(i) consist of indices of. outer R-splines. 
They correspond to complementJary index sets L(j) of indices of inner. B-spLimcs) 
i.e., i belongs to if and only ifl j belongs to J(i). For a given outer index j the 
indlex set T(j) is an array , i.e., a quadratic ot cubical arrangement of (n+ l) d inner 
indioes which is characterized by. a minimum distance to the index j. For, a given 

l*0f outer index j and an inner index i in the index set I^j), leti p\ be the d-var-iate 
polynomial ofl degree n in each variable which has tihe value 1 all the point i and at 
all other points of the array the value 0. Then the coupling coefficient Cij is 
given as the value of pi at the point j], i.e., = Pi^j). Thle specific values of the 
coupling aoofficianft can either, be tabulated for different degrees and relative 

15 positions of j and or can be easily computed using Lagrange polynomials. 
Figure 11 shows thle input and output data for this process step. 



Tf natural boundary conditions are given, the extended splines defined in equation 
(1) arc used without modification for further implementation of the FE process. Ml, 
20 on the oChcr hand, essential boundary) conditions are given, a weighting according 
t)o tlhc rule 

DM^^Biixl iS I (2) 

- has to. be performed. The pertinent interrogation takes place in an optional prooess 
step 6a. The functions defined in tihis- way are called weighted extended B-splines 
25 (jWEB-spIincs). Formally, the extended B-splines uaed under natural boundary 
conditions correspond to the spcoJaJ case w(x) = 1. They, are, therefore, also called 
WEB-spIincs. For the cage of esscntilaf boundary conditions, tibe weight function ?/; 
ia characterized as follows: Fbm all points x of the simulation region, vi(x] can be 
bounded from above and below, by positive! constants, wihich are independent ofl x, 
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tiimes the distance dist(x) of the point x from the boundary P. Tn other words, 10 
is positive within H and tends to aero in the vicinity of the boundary V as fasti as 
tine distance function dist. For simulation! regions wiltioh arc bounded hy 
elementary, geometrical objects (circles, planes, ellipses, etc.) a suitable weight 
5 function can optionally be given in explicit analytic formi. Otherwise, commutation, 
rules, should be used which 1 typically represent a smoothing of; the distance 
function. The scaling factor 1/wfoi) is calculated by evaluating the weight 
function at the weight point a^. Tlhis can be any point in the support of the 
B-spline ^ which its at least half th'e bound s/2 away from the boundary). 

lid 

As a result of the process of the present invention, a computation rule for the 
WEB-splines Bi (compare definitions (1) and (2)'j is obtained which have all 
favorable properties according to the objectives of thle invention. Thus the FIE 
method can he continued according to the prior art. Hut, in doinig so, it is possible 
15 to advantageously exploit the regular grid struoture of the basis functions 
according to claim 11. 

In summary, far the prooess of thle present invention the ooupling of outer and 
inner D-spline is, among othler things, important. As a result, the constructed 

20 basis has the properties wihioh are essential for FE computations. In particular, a 
basis Bi (ji from the index set /) according to the present invention, is stable, 
uniformly wrfth' respect to tho grid width h, and the error has the same order as for 
the R-spIincs &n ?» the approximation of smooth! functions which satisfy the same 
boundary conditions. On the other, hand, the fulfillment of essential boundary 

25 conditions is ensured by using the weJght function w. 

Other objects, advantages amid sailrerct features ofi the present invention willJ become 
apparent from flic following detailed description wihich, taken in conjunction with 
the drawings, listed below, disc-loses preferred embodiments of: the present 
3d invention. TIlc features mentioned in the claims and in the specification can We 
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essential to th'c invention individually or in any combination. 
Brief Description of the DrawSngs 

Referring to the drawings which form a paTt of this disclosure) 

Figure 1 is .a flow chlart of the individual steps of the prior, art in the course 

of the finite element simulation and integrates the determtinaitkm of 

the WEB basis into this prooess; 
Figure 2a compares certain finite elements of the prior art to the WEB 

element shown in Figure 2b and lists tlhe parameters relevant to 

finite element approximations: 
Figure 3 is a flow chart of the prooess steps for determining the W1EB basis; 
Figur.es 4a and 4b show a support and the corresponding tensor product 

B-spllne of degree 2; 
Figures 5a amd 5b illustrates the problem formulation of! the first 

embodiment (displacement of: a membrane under constant 

pressure) and of the corresponding solution; 
Figure 6 shows the ceil types for the first embodiment; 
Figure 7. sunveys the input and output data of the prooess for determining 

tlic cell typesj 

Figure 8 illustnartfes, by way of) example, the classification of the B-splines 

for tlic firwti embodiment]; 
Figure 9 outlines the input and output data of the prooess for classifying 

the D-aplines; 

Figure 10 shows the onuptftig coefficients of an outer B-spline and the 
corresponding inner B-splincs for the first embodiment; 

Figure 11 surveys the input and output data of the process for computing 
the coupling coefficient's) 



Figures 12a and' 12b illustrates tUc construction of tUc wciglilt function of 

the preferred embodiment; 
Figure 13 shows the support of: a. WER-spline and the corresponding 

coupling coefficients for the first embodiment] 
Figures 14a and 14b explains the problem formulation of: a second 

embodiment (incompressible flow) and its solution using the fiowi 

lines and th'e distribution of the flow velocity); 
Figtm&s 15a to 15c show the ohamging classification of B-spKnes for the 

B-spline degrees n = l 1 , 2, 3 and the same grid width for. the second 

embodiment) 

Figures 16a to U6c provide information about the error development in th'e 
finite element approximation using WEB-spliraes and about the 
computing time behavior of WEB approximations for the second 
embodiment) 

Figures 17a and U7b compare the WiEB basis to a prooess based on linear 

trial functions on a triangulation (priori art); and 
Figure 18 shows a computer system according to the present invention. 



- 11 - 



Detailed Description of the Invention 

One especially favorabl'e embodiment of the pnoccss of th'c present invention), called" 
the WEB prooess, te dleterrnined by the following specifications. 

Q The bound .s is chosen such that the inner B-splines are characterized by; rcquininig 
thlat at least ane of the grid cells of their support lies completely in tha simulation 
region Q. Since to determine the relevant B-splines, the intersection of: the gridl 
cellls and the boundary) T must be computed anyway, the classification requires rro 
significant additional computing time. The weight point) x t is chosen as the 
10 midpoint of a grid cell in the support of the B-spline hi whfch lies completely in 
tlli'n simulatilon region Q. This is also efficiently possible since the determination of 
one such cell is already, part of the classification routine. 

Tf no explicit analytic form of the weight function is knowin, it is defined by) 



Figures 12a and 12b illustrate the construction of the weight function. Here, the 
parameter & indicates the width of the strip Qg w?tihin which the weight function 
varies between the value 0 ot tihle boundary) ofl lihe simulation region and the value 
1 on the plateau am Q\Qs. The parameter <8 is chosen such that the smoothness of! 
20 the weight function is ensured. 

One important advantage of the prooess is that no meshing of the simulation 
region is necessary. Tn teahninal applications, this results in clear savings of 
computing time and storage capacity and simplifies the course of thle simulation. 
25 The process structure for two- and thrac-dimensianal problems is formally and 
tcchmcally largely identical. This enables time- and cost-saving implementations of 
solutions for cJivcrsc applications based on uniform prognarn structures. The use of 
Bl-splincs corresponds to f ho induisUrial standard in the modeling of geometrical 
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objects, and thus forms a natural connection between FE and GATING AM 
applications. Extensive existing program libraries from both fields can tic used for 
implementing a FE simulation based cm thle process of tlie present invention. The 
basis funotions constructed using the WEB process have all standard properties of 
5 finite elements. This includes especially the stability of the basis. T:t implies, for 
example, that for. linear elliptic boundary value problems the ooniditiom niumbor of 
the resulting system of equations does not grow faster than for optimal 
triangulations as the grid width becomes smaller. For. applications, this means, for 
example, that linear systems of! equations as they typically arise in FE methods, 
10 cam be efficiently solved by iterative algorithms. Furthermore, for a given degree, 
tlh'e approximation order is maximal and the number of necessary parameters 
minimal. Thus, very accurate approximations are possible with a relatively small 
mrmher of parameters. Specifically, this can mean that the accuracies which so far 
required the use of mainframe computers can now be achieved witih workstations. 
15 The regular grid structure of the basis of the present invention permits a very 
efficient implementation, especially for assembling and solving.FE systems. 
Moreover, by using the weight function, the boundary, conditions can be satisfied 
during simulation without affecting the regular grid structure of tfhe basis 
functions. Fmailly,. by using mult igrid methods to solve the linear systems arising 
20 iff linear elliptic boundary value problems, one oan achieve that the overall solution 
dime is proportional to the number of unknown coefficients, and thus optimal. 

Tlic process of the present invention in the special preferred embodiment (WTEB 
process), is illustrated using the firsti embodiment shown in Figures 5a and 9b. The 
25 differential equation and boundary conditions are chosen to be very elementary so 
tliat in addition to the construction of the WUB basis of the present invention the 
entire course of the FE simulation can he followed without major additional effort!. 

Figure 5a shows an clastic mambraTie : which is fixed along the edge F of a planar 
30 simulation region H, and on whinh a constant prassune / = 1; acts iaside the region. 



Wiflh suitable normialfalian, the displacement u satisfies tfhc Poisson equation 
w.LthJ homogeneous boundary] conditions, 

— Au = 1 in; H 
u — 0 on Fl. 

The displacement u or the deflection of the membrane is depicted in Figure 5b. As 
described above, the WIBB process is divided into the following steps. 

Input 1 of the simulation region O] The boundary F is a periodic spline curve of: 
liQ degree 6, which is stored by] its control points 20 (in Figure 5a identified with black' 

dots). 

Input 2 of the boundary conditions] the homogeneous boundary condition is 
essential so that the construction of a weight funotion is necessary. 

15 

Tnput 3 of. tho control parameters: The degree n = 2, and, in order to make th'e 
figures caster, to understand, a relatively large grid width h = U/3 are used. 

Dcticnninatlioii 4 of the cell types: As illustrated ifn Figure 6, the simulation region 
20 is covered by a grid 21. which contains the grid cells of the supports of all 

D-spKncs of potential relevance fbr thle basis construction. The type determination 
m the example yields 69 outer grid colls 22, and 11 inner grid cells 24, and 20 grid 
cells 23 on th'c boundary. 

Classification 3 of the B-splinesi Here, the support of the R-spline bf. is a square 
Qk with corners 

C*"if, + 3, k 2 )h, (kji +. 3, k 2 -H 3)/i, (h u + 3)rt; 

25 Q(-i,o) and Q(2,ii) arc shown in Figure ff. The grid points kh of the relevant 

Hl-splincs, for mlncli! Qn intersects tha ilntleriar of the simulation region, are marked 
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in Figure 8 by a point or a circle. Alt grid 1 points ih for inner R-spfincs $i from the 
index list I), for. which' at least one cell of. the support Qi lies entirely within ft, are 
marked by. a point*. For. example, for i = (—4. 0). th'e grid cell (—2, 0)h +: [0, h] 2 lies 
entirely in Q. All grid points jh for outer R-splines (j from the imdax list J), for 
5 which no cell of the support Qj lies entirely in f2, are marked by a circle. 

Computation 6 of the coupling coefficients: To determine the coupling coefficients 
eij fbr each fixed j of thle fndex list J the nearest 3 x 3-array 

= {ft, It + 1, 4 + 2J x {4, ft -H 1, 4 + 2} 

liO of indices tn f ; is sought. In Figure 10, for th'e outer grid point j — (-1!, 2), which is 
markled with a circle, th'e array is identified with points. Figure 10 likewise 
shows tihe corresponding coupling coefficients in a matrix representation. They are 
computed hiy bivariatie interpolation. Far example, the interpolating polynomial: 
fox i = (-1. —1) is 

15 Pi(x) = -xv(x, -H 2)x 2 (x 2 - l)/2. 

Its value at the point x — j — (— 1, 2) is p»(i) = 1 = e^. One notices thiat many, of 
tiUc coupling coefficients are 0. TWis is a typical phenomenon. The coupling 
coefficients Qij are not; equal; to 0 for all i of the index list Ifo) only- if the indices i 
arc different from tho 'index j. in each component. 

20 

Computation rule 7 for the weight function) The weight function is grvien by, 
equation J3] with n = 2 and S = 0.2. The parameter <S is computed numerically. It 
must be small enough' to prevent; singularities of the distance function. To compute 
tlVc distance function, the process generates a program which uses Newton's 
25 method. Since flic weight function is not eojual to 1 only in a boundary strip, the 
complexity in the subsequent evaluations is low. 

Output: Figure 13 shows the support of a WFB-spline B { and the data necessary 
for its description. These arc the index list; J§i) orfl trie outer R-splines bj coupled 
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with' ft, the coupling coefficients Cij, and the weight' paint x t . These data arc used 
in conjunction with th'e weight function for- generating the computation rule for 
the WEB-splines. 

5 The further- course of tihe FE simulation follows the prior art. 

Assembly: 9 of: tihe FE system: The entries of the system matrix and of the 
right-hand 1 side are 



J a Jn 

ffl The system of equations GQ = F for th'e basis coefficients C t - in this example has 
dimension: 31. The matrix entries dj are computed using numerical integration, 
l?kcwtfge the integrals F fc . 

Solution 10 of the FE system: The Galerkin system is solved iteratively with the 
15 conjugate gradient method with SSOR preconditioning used to accelerate 
convergence. After 24 iteration steps the solution is found within machine 
accuracy (jtofer-anoe <] le — 14). 

Computation 11 fmd: outlput 12 of. the approximation] The approximation 
20 computed with the process as claimed in th'e invention 3s u — Yli QB% an< * ts 
shown graphically ini Figure 5b. The relative error of th!e Z^-norm is 0.028. 

The efficiency of the process of the present invention in the special preferred 
embodiment (WEB process) is rllustratied in a second embodiment! using tihe 
25 simulation of an incompressible flow. The arrangement of two circular obstacles 
shawm in Flfgurc 14a in a channel wiflth parallel boundaries serves to illustrate the 
principal strategy. Flor complicated geometries, as are typically present in specific 
applications, the process workis completely analogously and efficiently. Tn 
Fligurc 14a the stream lines 25 ame shown within the region bounded by Ti, to r 4 




k, i € A 



- 16 - 



and by P 5 andl F! 6 . The differential equation is: 



An - 0 



In O 



with th'e boundary conditions 



— — = vo on 



u i 
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5 Thle flow velocity v == — grad u is shown in the bottom half of thle figure. 

The construction of the WEB basis of the present invention proceeds completely, 
analogously; to the first embodiment. The solo difference is that a weight function 
is not necessary because of the natural boundary! conditions. 



Figures l'oa to U5c show the classification of the relevant B-splines for different 
degrees n (sec also Figure 8). In the figure, the inner B-splines which are taken 
into the WRR basils without extension, are marked by solid triangles. For small h 
IJhc number of these B-splines increases, i.e., Di = ft for most of! the WEB basis. 
15 For degree n = 3 this is the case for 236 of 252 indices in the example. 

Figure 16a shows in two diagrams the numerically determined relative X^-error ofl 
the potential (left half of: the figure) as a function of the grid width H = 2~* with 
k =« 1, . . . , 5 and the numerically) estimated order of convergence m (bright halfi ofl 

2Q the figure). Here, for different degrees of the WEB-spIine the fallowing markers aTe 
used: * (!jt 1), o (n = 2). a = 3)', a (n = 4) and « (n = 5j. As expected, 
m « it +' 1, i.e., air approximate error reduction by a factor 2 n+l when the grid 
width is cut in half. Analogously, for the relative approximation error of; the flow 
velocity shown in Figure l'6b (i/^-narm of the solution, left half ofl the figure), am 

25 order of convergence m « n {right balf of figure) is obtained with an associated 
error reduction bfy roughly ai factor 2" wihen thle grid width is cut in halfl 



10 



Figure 16c (right half of the figure) shows the computing time in seconds for 
construction! of the WEB basis as a function af the number of resulting basis 
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functions, measured on a Ftenftium TT processor witl»l 400 MHIz. For example, for 
construction of a WEB basis of degree 3 with' grrd width h = 0.12ft with 2726 
WEB-splines 1.32 seconds are necessary. One notices that tli'e complexity far 
generating th'e WEB basis is largely independent of the degree n of the basis. Tn 
5 flhe left half of Figure 16c the number of CG-itoraitions relative to r,bo number of 
basis functions is shown. Thus, for tihle corresponding system with 2726 unknowns, 
65 POG-iterations are required. The total computing time including assembling 
and solving the Qalerkin system is roughly 2.48 seconds. 

10 Figures 17a and 17b compare tWe WEB process with a standard solution prooess 
wWicb meshes or triangulates the simulation region ^Fig. U7a) and uses hat 
functions. The graph shows in Fig. 17b the L 2 -error relative to the number of 
parameters. The results of. the standard solver are marked with boldfaced 
diamonds and are compared to the results achieved using tJhe WEB basis of 

15 degrees U tlo 5. Fbr example, an accuracy of 10" 2 is achfieved with the WEB 

process by using 213 basis functions witih degree 2 and an overall computer time of 
0.6 seconds. To achieve the same accuracy^ the standard method with) linear hat 
functions rcquirod A657 basis functions. 

2Q In the assessment, of the standard solution process two other aspects must be 
considered. On the otlfer hand, Figure 1 71b illustrates that even a moderate 
accuracy of 10~ 3 cam only He achieved with hat functions when far more than one 
million coefficients are used. This shbws that wiWen using hat functions accurate 
results generally rcqurne an enormous competing and storage capacity or. cannot be 

25 achieved at all with the prior art. On tWe other hand, the complexity required for 
mcsliing increases with" the complexity of the simulation region. Tn contrast to 
realistic applications, flic region studied Were is comparatively simple to 
triangulate due to. it's simple stirunture. 

30 The two-dimensional example shows the performance gain by the WEB prooess. 
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An even greatcn increase in performance is possible in three-dimensional problems. 
On the one Wand, the complexity for rmeshing, which' is eliminated: in the WEB 
process, is much greater. On the other hand, the reduction in the number of 
required basis functions becomes rnuoh more noticeable than in the 
5 two-dimensional case. 

Flilgure 18 shows a device according to the present invention, especially a computer 
system 30, with input devices 31, 32, 33, output devices 34 and a control unit 35 
which controls the course of tJhe process. To carry out the prooess of the present 

10 invention and in particular for purposes of! paraJlelizatkm of the pertinent 

computations, the central control unit 35 preferably uses several arithmetic logic 
units (|ALU) or even several central processing units (CPU) 36. These allow 
especially parallel processing for, the process steps classification 5* of the B-splines, 
in partiauliar also intersection of the regular grid with the simulation region fi, 

15 determination 6 of the coupling coefficients e^, and/or evaluation of the weighti 
function w(x) at points x of the simulation region Q. 

The computer, units 36 here access the common data resources of the storage unit 
371. The data aain be iniput, for example, by) a keyboard 31, a machine-readable 

20 data medium 38 viia a corresponding read station 32 and/|or via a wire or wireless 
data nctftvork wirth a receiver station 33. Via the read station 32 or a pertinent 
data medium 38, the control program, which controls the process execution, can be 
input, and, for example, can bo permanently filed on the storage media 37. 
Accordingly, th'c output devices 34 can be a printer, a monitOT, a write station for 

25 a> machine- readable data medium arcd/on the transmitting station of a wire or 
wireless data- network. 

While various embodiments have been chosen to illustrate the invention, it will be 
understood hy those skilled in tJhc nrt that various changes and modifications can 
30 he made (ihcrcin without departing from the scope of the invention as defined in 



